Measurement-Based Quantum Computation on Symmetry Breaking Thermal States 
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We consider measurement-based quantum computation (MBQC) on thermal states of the inter- 
acting cluster Hamiltonian containing interactions between the cluster stabilizers that undergoes 
thermal phase transitions. We show that the long-range order of the symmetry breaking thermal 
states below a critical temperature drastically enhance the robustness of MBQC against thermal 
excitations. Specifically, we show the enhancement in two-dimensional cases and prove that MBQC 
is topologically protected below the critical temperature in three-dimensional cases. The interact- 
ing cluster Hamiltonian allows us to perform MBQC even at a temperature an order of magnitude 
higher than that of the free cluster Hamiltonian. 



Introduction. — Measurement-based quantum compu- 
tation (MBQC) is a paradigm for quantum computation, 
where a many-body entangled state is prepared as a uni- 
versal resource state, and quantum computation can be 
executed by adaptive single-qubit measurements on it [l[ . 
This paradigm provides a good clue to understand re- 
quirements on a system as a resource for universal quan- 
tum computation, making a bridge between quantum in- 
formation science and many-body physics. A central is- 
sue in this approach is to specify a many-body system 
whose ground or low-temperature thermal state can serve 
as a universal resource for MBQC. 

Ground states of several Hamiltonians such as the clus- 
ter Hamiltonians 0-(3] and valence-bond solid Hamilto- 
nians 0-0, have been found to be universal. At finite 
temperature, thermal states of several Hamiltonians have 
been shown to be useful as universal resources by protect- 
ing quantum information from errors originating from 
the thermal excitation by using quantum error correction 
However, these Hamiltonians do not undergo 
any physical (thermal nor quantum) phase transitions, 
although they exhibit a transition in computational ca- 
pability of MBQC by varying temperature. Thus we ad- 
dress a question whether it is possible to enhance the ro- 
bustness of MBQC against thermal excitations by intro- 
ducing a Hamiltonian that undergoes a phase transition. 
This can strengthen the connection between quantum in- 
formation science and many-body physics and can pro- 
vide an approach to understand the robustness of MBQC 
in terms of many-body physics. 

In this letter, we show that robustness of MBQC on 
thermal states can be enhanced drastically by introduc- 
ing interactions between the cluster stabilizers. The 
proposed Hamiltonian -interacting cluster Hamiltonian 
(iCH)- is transformed into a ferromagnetic Ising Hamil- 
tonian by unitary transformations. Hence, it undergoes 
a phase transition on two or higher dimensional lattices, 
leading to the symmetry breaking of the thermal states. 
By virtue of the ferromagnetic-type long-range order of 
such symmetry breaking states, MBQC becomes robust 



below the critical temperature. We first demonstrate 
this on a two-dimensional (2D) lattice and show that 
the fidelity of MBQC can be drastically improved below 
the critical temperature due to the long-range order, al- 
though it is not sufficiently large at the temperature just 
below the critical temperature. We further investigate 
topologicallyprotected MBQC on a three-dimensional 
(3D) lattice 0, Il3l - fl9| in order to achieve quantum com- 
putation of arbitrary accuracy at any temperature below 
the critical temperature. We show that the threshold 
value for the topologically protected MBQC is exactly 
equal to the critical temperature. Compared to the pre- 
vious Hamiltonian without the interactions between the 
cluster stabilizers 0, the temperature required for topo- 
logically protected MBQC is relaxed by more than one 
order of magnitude. 

Cluster Hamiltonian. — The cluster stabilizer on a lat- 
tice T is given by Ki = Xi ®j^y. Zj for each site i, where 
Ai (A = X, Y, Z) are the Pauli operators on the zth qubit, 
and Vi denotes the set of the vertices that are adjacent 
to the site i in the lattice T [J. The cluster state on 
the lattice T, |^r)i is defined by the simultaneous eigen- 
state of all cluster stabilizers with eigenvalue +1. The 
cluster Hamiltonian is defined by using the cluster stabi- 
lizers as Hi c = —JJ2i Ki ( see Fig- DI a ))> where J is 
a coupling constant. It is obvious that this Hamiltonian 
has the cluster state \^t) as its ground state. By using 
a unitary transformation Uf, the products of controllcd- 
Z gates on all bonds of the lattice T, the Hamiltonian 
can be transformed into an interaction-free Hamiltonian 
U-rHf c U-T = —J^iXi. We call this Hamiltonian a free 
cluster Hamiltonian (fCH) hereafter. It should be em- 
phasized that this system does not exhibit any phase 
transition since the thermal state is equivalent to that 
for the interaction-free Hamiltonian. 

The thermal state can be calculated as (FT Si) |^ r )(^ r |, 
where £ iP = (p+ e~ 2 ^ J Z iP Zi)/(l + e~ 2 ^ J ) f|, which can 
be regarded as a cluster state with an independent Z 
error on each qubit with probability ppj = e~ 2l3J /(l + 
e~ 2 P J ). It has been known that it is possible to per- 



2 




FIG. 1: (a) The fCH on the square lattice, (b) The iCH on 
the square lattice, (b) The interacting cluster Hamiltonian 
on the RHG lattice, (c) The 3D cluster state on the cubic 
lattice, (d) The gray-shaded qubits are measured in the Z 
basis to obtain the cluster state on the RHG lattice. 



form MBQC in a topologically protected way with a 3D 
cluster state on the so-called Rausscndorf-Harrington- 
Goyal (RHG) lattice [TMl- The threshold values have 
been numerically calculated to be ppj = 2.9-3.3% (i.e. 
T c = l/((3 c J)= 0.57-0.59) depending on the decoding 
algorithms [l5|, EH ■ These threshold values correspond 
to the critical concentrations of antifcrromagnetic inter- 
action at zero temperature and at the multicritical point 
(on the Nishimori line 2C |]_of the random plaqucttc Z2 
gauge model (RPGM) 



Next we introduce an interacting cluster Hamiltonian 
(iCH), Hi c = —jy)/jj\ KiKj where the summation runs 
over all bonds (ij) of the lattice T, and thus each cluster 
stabilizer interacts with its nearest neighbors as shown 
in Fig. [1] (b) for the case of the square lattice. This in- 
dicates that the iCH generally contains mutually depen- 
dent stabilizer operators. The iCH can be transformed 
to the ferromagnetic Ising Hamiltonian on the lattice T 
by Uj- as U-rH ic Uj- = —JJ2{ij) XiXj, which is denoted 
by i/ising- If the geometrical structure of the lattice is 
chosen properly, for example, two or higher dimensional 
lattices, phase transitions happen. Although each eigen- 
state of the iCH, |*), is degenerate with Ilier ^l^) ^ ue 
to the symmetry of the iCH, the symmetry breaking takes 
place below the critical temperature. Hence, the popu- 
lation of either |\&) or ILeT^! 1 ^) becomes much larger 
than that of the another, which results in a symmetry 
breaking thermal state exhibiting a long-range order. We 
can identify which of the symmetry breaking states is re- 
alized by measuring a cluster stabilizer. MBQC on such 
symmetry breaking thermal states would exploit robust- 
ness of long-range order. 

MBQC on the square lattice iCH. — Let us first con- 
sider MBQC in the iCH on a square lattice where a pe- 
riodic boundary condition is assumed. Since the iCH is 
unitarily equivalent to the Ising Hamiltonian, the sys- 
tem undergoes a phase transition at the critical temper- 
ature T c /J = 2/ln[l + s/2] [Uj]. To check whether the 
symmetry breaking thermal state leads to the robust- 
ness of MBQC, we consider to perform the identity and 



Hadamard gates. 

For performing the identity and Hadamard gates, we 
remove unnecessary qubits on the square lattice by mea- 
suring them in the Z basis and obtain a cluster state on 
a line pj. To transfer a state of the first qubit on the 
line to the Z-th qubit, we measure all qubits up to the 
(I — l)-th qubit in the X basis and obtain the measure- 
ment outcomes m = {m^} <G {0, l}*^ -1 ). Then we apply 
the Pauli bi-product to the Z-th qubit, B m = X r ™Z r ™ 

where r* = Yh=i m 2i and rf n = Y^ilx "*2*-i- The 
summation is taken by modulo 2 and \a~\ is the smallest 
integer larger or equal to a. When the number of the 
qubits that are measured in the X basis is even (odd), 
the identity (Hadamard) gate is implemented 

When the MBQC is performed on the thermal states, 
each measurement outcome includes an error originated 
from thermal excitations. By denoting such an error by 
e = {e,} £ {0, 1} X " -1 ), the z-th measurement outcome 
containing the error is given by m, = m, + ej (mod 2) 
where mi is the measurement outcome for the ideal clus- 
ter state. This results in an error B e of the gate operation 
through the Pauli bi-product B,~ n = B e B m . When the 
error is correlated so that both of r£ and rf are even, 
B e = I, that is, the gate operation does not include any 
errors. Thus, the probability of implementing the gates 
without any errors is equal to the overlap between the 
thermal state and the projector onto the subspace where 
both of r* and the are even. Such a probability is 
referred to as a gate fidelity and is given by 



F(l) = Tr 



n 



r(j-i)/2i 



2? 



-Pth 



(1) 



where pth is a thermal state of a given Hamiltonian H , 
pth = e~P H /Tre~@ H . The gate fidelity takes values be- 
tween 1/4 and 1, and the minimum gate fidelity 1/4 im- 
plies that the output state is the completely mixed state 
and the gate operations fail. In the case of fCH and iCH, 
by applying the unitary transformation Uj-, the gate fi- 
delity F(l) is expressed in terms of the many-body cor- 
relation functions of the Ising model. The results for 
the identity gates are shown in Fig. [5] for fCH and iCH 
(for the detailed calculations and the gate fidelities of the 
Hadamard gate derived analytically, see Appendix A and 
B). 

For fCH, the gate fidelities exponentially decrease with 
increase of temperature to 1/4 for any distance I. For 
iCH, the gate fidelities change differently below/above 
the critical temperature, namely, the gate fidelities also 
exhibit a transition at the critical temperature. Above 
the critical temperature, the gate fidelities arc close to 
1/4. In contrast, the symmetry braking thermal state 
appearing below the critical temperature leads dramatic 
improvement of the gate fidelities by exploiting its long- 
range order. The temperature required to perform the 
gate operation reliably for iCH is much higher than that 
for fCH, e.g., the gate fidelity for I = 8 is almost 1 if 
T < 1.0J for iCH and T < 0.3J for fCH. However, even 
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FIG. 2: The gate fidelities of the identity gates for various 
distance /. Fig. (a) shows the gate fidelity for fCH, F(2l) 
for I = 2 (x), / = 4 (□), I = 6 (+) and I = 8 (o). Fig. (b) 
shows the gate fidelity for iCH, F(2l) for I = 2 (x), I = 4 (□), 
i = 6 (+) and i = 8 (o). The vertical dashed line shows the 
critical temperature for the 2D iCH, T c / J — 2/ ln[l + v2] and 
the horizontal dashed line shows the minimum gate fidelity 
F = 1/4. 



for iCH, the fidelities just below the critical tempera- 
ture is not large enough to reliably carry out the gate 
operations, which motivates us to consider the topolog- 
ically protected MBQC on the symmetry breaking ther- 
mal states on a 3D lattice. 

Topologically protected MBQC. — Topologically pro- 
tected MBQC can be p erformed with the cluster state 
on the RHG lattice [lM^. The RHG lattice is denned 
by the set of cubes Q, the set of faces F q on each cube 
q € Q, and the set of edges Ef on each face f G F q . The 
qubits are located on each face and edge constituting the 
3D cluster state, where the nearest-neighbor stabilizer 
generators interact with each other, as shown in Fig. Q] 
(c). Similarly to the original case (17T - ll9| . the 3D cluster- 
state is subdivided into three regions, vacuum V, defect 
D, and singular qubits S used for the topological protec- 
tion, performing Clifford gates, and performing arbitrary 
single qubit gate, respectively. In the following, we con- 
sider the threshold value for the topological protection in 
the vacuum region V, since it solely determines a thresh- 
old of quantum computation [rjj - [l9| . 

In the vacuum region V, all qubits are measured in 
the X basis. Since II/eF = 11/gf tne P arr ty of 
the measurement outcomes on the six face qubits on each 
unit cell has to be even for the ideal cluster state. Thus 
the Z errors on the face qubits, say error chain C, are 
detected at the boundary dC, which is called the error 
syndrome of C, since the error chain C anticommutcs 
with Il/eF there. This is also the case for the Z 
errors on the edge qubits, say error chain C, since the 
edge qubits are the face qubits on the dual lattice. In the 
case of the fCH on the RHG lattice 0, the error chains 
C and C, denoted by C = (C, C) are not correlated and 
can be treated independently. However, in the case of the 
iCH, the primal and dual error chains of C are strongly 
correlated and have to be treated simultaneously. 

By using the error syndrome dC, we infer the actual 
locations of the errors. To this end, the probability of 



a hypothetical error chain C, which has the same error 
syndrome dC, is calculated to be 



P (C'\dC) = Af~ 



exp 



ufuf 



<//> 



dC'=dC 



(2) 



where AT is the normalization factor. We used the knowl- 
edge that the errors occur with a ferromagnetic Ising-type 
distribution, which is characterized by a parameter /?' in- 
dependently of the physical inverse temperature ft. The 

-1 for / e C 



indicator function is defined as 



and 



1 for / ^ C", specifying the location of the er- 
rors. Since dC — dC, we have C = C+C for trivial loops 
(cycles) C = (L,L), where L is a trivial loop for the lat- 
tice and L is for the dual lattice, such that dC = 0. In or- 
der to solve the loop condition, we introduce gauge vari- 
ables on the edges of primal and dual lattices defined by 



P 



I 



UeeE f a e and P f 



n 



In this 



parameterization, II/eF Pf = 1 an< ^ ®L — is automati- 
cally satisfied. As a result we obtain the Gibbs-Boltzman 
distribution p(C'\dC) = ^f-i e -0'Hc(<r,s) under a Hamil- 

tonian given by Hc{(J, <r) = — J 



,c„,c 



jPfPt, which 



{f}) u f u jnn 
we call correlated RPGM (cRPGM). The sign of the two 
plaquette interaction UjU^ representing the randomness 
of the model is determined by the actual error chain C 

with the distribution p(C) = W~ 1 e' 3J5: <//> param- 
eterized by the physical inverse temperature (3. 

Since the threshold value for topologically protected 
MBQC corresponds to the critical point of the cRPGM 
[3], our goal is to identify it. Let us consider the optimal 
case of f3' = where the actual and hypothetical error 
distributions are the same. This condition is referred to 
as the Nishimori line poj in spin glass theory. In this 
case, the internal energy is given by 



[(H c (a,a)) th ]c = ^>(C) £ 



ZcW) 



}th denotes the thermal average and Zc{(3) = 
-PHc(cr,a) - 1S j-jjg p ar tition function. We take 



where 

{(T e ,0 

the ensemble average of the error distributions [-]c since 
the cRPGM has self-averaging property [22|. With the 
aid of the gauge symmetry [20( , the Hamiltonian He (<r, a) 
is invariant under the following gauge transformations, 



where P'j 



,Cpl 



a e -> er e cr e , uj -> ujfp erg ~> o- g a' s , 

UeeE f a 'e and Pf = riee% ^S- 0n tllC ° ther 

hand, these transformations changes the distribution 
p(C) into Af~ 1 eP J ^<f-f> u f u ? p 'f p 'f = -i e -f>H e (*',s f ) = 
p'(C), which corresponds to the Gibbs-Boltzman distri- 
bution for the cRPGM. Since £ {(T ,. ff , } p'(C) = Z c ((3), we 
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obtain 

= AT 1 ]Ttfi sing e-' 3 ^»- = (H lsing ) th , 
c 

where \C\ is the number of the loop configurations, 
-fusing = —JJ2(f J) u^u^, and A/" is defined in Eq. J2]) as 
the partition function of the Ising model. For the trans- 
formation from the third to forth lines, we take the sum- 
mation ~^2 a - = \C\ by using the fact that u'fPf = u / +i 

and u< fPf = U< f +L with trivial loops C and ^ c = J2c+c- 
Thus, the internal energy of cRPGM is equivalent to that 
of the Ising model without any randomness. 

In the Ising model on the RHG lattice, the internal 
energy has a non-analytical point at T c — 2.8, which is 
evaluated by the exchange Monte Carlo simulation (23| . 
Therefore we can conclude that the internal energy of 
the cRPGM along the Nishimori line also has a non- 
analytical point at T c = 2.8, which is the phase bound- 
ary of the Higgs (ordered) and confinement (disordered) 
phases [14T - [l6| . In the Higgs phase, the loop configura- 
tions C of large perimeters are exponentially suppressed. 
Thus the logical error probability, which is characterized 
by the loop configurations of non-trivial topology, is de- 
creased exponentially by increasing the size of the sys- 
tem (see Appendix C for the decoding methods). That 
is, the transition point of the performance in topologi- 
cally protected MBQC on the symmetry breaking ther- 
mal states exactly determined by the critical temperature 
of the phase transition in the underlying physical system. 

The cluster state on the RHG lattice can be also ob- 
tained from other lattices such as simple cubic (SC), face- 
centered cubic (FCC), close-pack hexagonal (CPH) lat- 
tices by measuring appropriate qubits in the Z basis as 
shown in Fig. [TJ (d) in case for the SC lattice. Since the 
thermal errors commute with the Z basis measurements 
in our model, they do not induce any additional errors. 
Also in these cases, by using the gauge transformation, 
we can show that the thresholds for topologically pro- 
tected MBQC again given by the critical temperatures 
of the iCHs on those lattices. The critical temperatures 
of the Ising models on the SC, FCC, and CPH lattices 
have been calculated numerically as T c — 4.5, 9.3, and 
9.8, respectively (24[, which arc higher than T c = 2.8 
for the RHG lattice, since each site interacts with more 
neighboring sites. In comparison with T c = 0.59 for the 
fCH, the iCHs with the long-range order relaxes the tem- 
perature required for topologically protected MBQC by 
more than one order of magnitude. In the fCHs, the lat- 
tice structures do not change the threshold value for the 
topological protection since the thermal errors occur in- 



dependently for each qubit. Contrary in the iCHs, the 
underlying lattice structures take a very important role 
in robustness against the thermal excitation by making 
use of physical cooperative phenomena. 

Conclusions and Discussions. — We have demon- 
strated that physical cooperative phenomena of a system 
can help MBQC on the system even at finite tempera- 
ture. We have first shown that, in a square lattice, the 
gate fidelities of the identity gates for iCH are drastically 
improved compared to those for fCH below the critical 
temperature. It has been also shown that the fidelities 
are not sufficiently large for performing MBQC reliably at 
the temperature just below the critical temperature even 
for iCH. In the 3D cases, MBQC on the thermal states are 
topologically protected below the critical temperatures of 
the underlying physical system, which allows us to per- 
form MBQC on the symmetry breaking thermal states 
even at much higher temperatures than the models with- 
out physical cooperative phenomena. A promising way 
to design these many-body interactions used in both fCH 
and iCH is the stabilizer pumping scheme [251428T ] . Al- 
though achieving larger many-body interactions requires 
more unitary operations in the scheme, the required tem- 
peratures for performing topologically protected MBQC 
is significantly relaxed for iCH. 

In the present work, we have considered only the Ising- 
type interaction in the stabilizer Hamiltonian. We can 
also construct the iCHs, which are unitarily equivalent 
to other spin models such as the Potts, XY, and the 
Heisenberg models. It is an interesting future work to 
study the relation between the ordered phase and quan- 
tum information tasks in such models. This will open 
up a new approach to make use of physical cooperative 
phenomena for quantum information processing. 
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Appendix A: Correlation functions in the 
interacting cluster Hamiltonians 

We show the detailed calculation of correlation func- 
tions in the interacting cluster Hamiltonians (iCHs) on a 
square lattice. The Hamiltonian is given by 

H ic = -jJ2KiKj, (Al) 

<«) 

where the summation runs over all nearest neighbor 
bonds, (ij), of the square lattice T. The thermal state 
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of the iCH at temperature T is denoted by p t h = 
e~P Hic / Z ic , where [3 := 1/T is the inverse temperature, 
and Zi c = Ti-e~P Hic is the partition function. 

The fc-body correlation functions at temperature T is 
defined by 

k k 

(II *(W»)> : = Tr II K (**Jn) ■ ■ ■ K (i k ,3 k )P^, (A2) 
n— 1 n— 1 

where (i n ,j n ) denotes a coordinate of the square lattice. 
For simplicity, we assume i n < i m and j n < j m for n < 

TO. 

By the unitary transformation Uf, which is a product 
of controlled-Z gates on all bonds, iCH is transformed 
into the same form as the Hamiltonian of the Ising model, 
UrH ic Uj- = -JJ2(i 3 ) XiXj = Rising- Hence, the corre- 
lation functions are expressed in terms of those of the 
Ising model: 



(A3) 



where (^4)i sing = TTAe~^ HlBine / Zi sing . In the following, 
we present the fc-body correlation functions in the square- 
lattice Ising Hamiltonian for even and odd k. 



1. Even-body correlation functions 

We present the exact formulas of even-body correla- 
tion functions on the i-th row, (H„ = i -X7i,j„))ismg- Since 
the Hamiltonian has a translational invariance, we con- 
sider only the correlation functions on the first row, 
i = l, without loss of generality and we simply denote 
the correlation functions by (FIn=i -^jn) ising- The cor- 
relation functions can be analytic ally obtained by using 
the method of a transfer matrix (291 ], and the resultant 
exact formulae are given as follows. 

We define a function 



c(6) 



2z(l + z 2 ) - z 2 (l - z 2 )e l0 - (1 - z 2 )e 



y/[(l + z 2 ) 2 - 2z(l - z 2 ) cos<9] 2 - 4z 2 (l - z 2 ) 2 ' 

(A4) 

where z :— tanh^J, and its Fourier transformation 



( r 



i 

2^ 



d0e- trO c(9) 



(A5) 



Using C r , the correlation functions in the thermodynamic 
limit are analytically given by 

2k k hn-l 

n—1 cr n—1 l—j^n — i 

(A6) 

where tr is an element of a permutation group defined by 

k 

a e Perm (J {j 2 „-i + l,j 2n -i + 2, • • • , j 2 „}. (A7) 




(b) 
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FIG. 3: The numerical results and the analytical results for 

(a) a correlation function (Xn^X^i^Xn^^ti.^isiag and 

(b) a correlation function (]~Ife=i ^(i,2fc-i))ism g - The lines 
represent the exact results and the points (o) represent the 
numerical results. In Fig. (b), the red vertical dashed line 
shows the critical temperature given by T c = 2/ln[l + x/2]- 



2. Odd-body correlation functions 

In contrast to the even-body correlation functions, the 
odd-body correlation functions are zero because of the 
symmetry of the Ising model. To see this, apply a flip op- 
erator to all qubits, namely, a product of all Z operators 
on each site. Since the flip operator commutes with the 
Ising Hamiltonian i?i S ing and anti-commutes with Xuj\ 
for any we obtain 



2k+l 



2k+l 



ing, 



leading to 



2k+l 



71=1 



,) /Ising 



(A8) 



(A9) 



Note that, when the symmetry of the system is broken, 
the odd-body correlation functions are non-zero. The 
odd-body correlation functions for symmetry breaking 
states have been analytically calculated under a certain 
assumption (30j . Here, we numerically calculate them by 
the Monte-Carlo simulation without using the assump- 
tion since it has not been rigorously shown whether the 
assumption always holds, up to our knowledge. Since 
the initial state of the Monte-Carlo simulation is break- 
ing the symmetry, it results in the correlation functions 
for the symmetry breaking states. 



3. Monte-Carlo simulation for correlation functions 

We provide numerical results of the odd-body correla- 
tion functions obtained by the Monte-Carlo simulation. 
The Monte-Carlo simulation performed with the lattice 
size 150 x 150 and the number of sampling 10 5 . To 
check the finite size effects of the simulation, we prelim- 
inarily compare the numerical results of the even-body 
correlation functions with the analytical results given in 
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FIG. 4: The numerical results and the analytical results for 
the odd-body correlation functions (IlntJi 1 -^(l,n))lsi»g for (a) 
k — 1 , (b) k = 2, (c) k = 3 and (d) k = 4. In each figure, 
the red vertical dashed line shows the critical temperature 
T c = 2/ ln[l + ^]. 
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Sec. IA 11 and confirm that the finite size effects are neg- 
ligible in this case (see Fig. [3]). 

The numerical results for the odd-body correlation 
functions, (Jlnta 1 -^(i.n))ising, arc shown in Fig. |4] The 
correlation functions are zero above the critical temper- 
ature. Below the critical temperature, the correlation 
functions are non-zero and are decreasing with increas- 
ing temperature. Note that, for larger k, the correlations 
decrease more quickly with increasing temperature, im- 
plying that the many-body correlations are not robust 
against the thermal excitations. The correlations at suf- 
ficiently low temperature are almost unity independently 
of k, that is, the spins are in an ordered direction. 



Appendix B: Gate fidelities of the Hadamard gate 

Based on the correlation functions, we calculate the 
gate fidelity of the Hadamard gate acting on the qubits 
at (1, 1) and (1, 21 + 1) where I = 1, 2, • • • . When a state 
is a thermal state pth, the gate fidelity is given by 



0.8 



F(2l + 1) = Tr 



I 



n 



= 1 K (l,2i) 



Plh 



For a thermal state of a free cluster Hamiltonian (fCH), 
denned by Hf c = —JJ2i Ki> ^ is straightforward to cal- 
culate F(2l + 1) since thermal excitations on each qubit 
occur independently: F fc (2l + 1) = (1 + tanh' (3 J) 2 /4- 
For a thermal state of iCH, the gate fidelity can be cal- 
culated by using the results in Sec.[XJ In particular, when 
I is even, the gate fidelity F(2l + 1) contains only even- 
body correlation functions so that it can be analytically 



1.6 2.4 
Temperature 



3.2 



FIG. 5: The gate fidelity of the Hadamard gates for various 
distance I. Fig. (a) shows the gate fidelity for fCH, F(2l + I) 
for I = 2 (x), I = 4 (□), I = 6 (+) and I = 8 (o). Fig. (b) 
shows the gate fidelity for iCH, F(2l + 1) for / = 2 (x), I = 4 
(□), I — 6 (+) and 1 = 8 (o). The vertical dashed line shows 
the critical temperature for the 2D iCH, T c /J = 2/ln[l + 
\/2] and the horizontal dashed line shows the minimum gate 
fidelity F = 1/4. 



obtained by using Eq. (|A6|) . The results are shown in 
Fig. El 

Compared to the exponential decay for fCH, the gate 
fidelity for iCH changes its behavior at the critical tem- 
perature. This transition of the gate fidelity is clearly 
observed by the derivative of the gate fidelity as shown in 
Fig. [6] Similarly to the gate fidelity of the identity gate, 
the gate fidelity of the Hadamard gate for iCH is drasti- 
cally improved below the critical temperature, although 
they are not sufficiently large for performing MBQC re- 
liably at the temperature just below the critical temper- 
ature. 



Appendix C: Decoding methods 

On the vacuum region V of the RHG lattice, we per- 
form X basis measurements for topological quantum er- 
ror correction The thermal excitation results in 
the Z errors on the face qubits on the primal and dual 
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FIG. 6: The first derivative of the gate fidelity of the 
Hadamard gates for I = 2 (x), I = 4 (□), I = 6 (+) and 
I = 8 (a). The vertical dashed line shows the critical temper- 
ature for the 2D iCH, T c /J = 2/ ln[l + \/2]. 



RHG lattices, which arc associated with the edges on 
the dual and primal lattices respectively, and denoted by 
a chain C = (C,C). The Z errors are detected at the 
boundary DC of the error chain C, since H/gf be- 
comes odd there. In the decoding procedure, we infer 
the location of the error chain C from the error syndrome 
dC. 

In the case of fCH on the RHG lattice, the Z errors 
due to the thermal excitations occur independently for 
each qubit. Thus the minimum- weight-perfect-matching 
algorithm (MWPMA) can be used as a suboptimal but 
good decoding method, resulting in the threshold value 
T c = 1/(P C J) = 0.57 (ppj = 2.9%) pi, which is close to 
the optimal one T c = l/(/3 c J) = 0.59 (jppj = 3.3%) [H]. 
In the case of iCH, the MWPM is far from optimal, since 
the Z errors are strongly correlated due to the interaction 
terms KiKj. This is numerically observed in Fig. [7] The 
Fig. [7] provides the logical error probability, with which 
the error correction is failed, as a function of temperature 
T = l/(j3 J). If the temperature is smaller than T 
1 .!). the logical error probability exponentially decreases 
with the lattice size TV, which implies that the decoding 
succeeds. However, this threshold T = 1.9 is smaller than 
the optimal temperature T c = 2.8 given by the critical 
temperature of the Ising model on the RHG lattice. 

To achieve a successful decoding at the temperature 
close to the optimal temperature, we should take into 
account the correlations between the Z errors. One 
way is making use of the free energy (3Fc+Vi(P) = 
— \nZc+Vi (/?) °f the cRPGM for hypothetical error 
chains C + Vi. Here C is an arbitrary error chain that 
satisfies dC = dC, and Vi is a logical operator, which is 



represented by a non-trivial loop belonging to a homol- 
ogy class i [lj]. In the decoding process, if the actual 
error chain C and estimated error chain C + Vi belong to 
the same homology class, the error correction succeeds. 
To succeed the error correction with a high probability, 



we infer the most likely homology class from the error 
syndrome dC. 

1 




1.5 1.6 1.7 1.8 1.9 2 2.1 2.2 

Temperature 

FIG. 7: The logical error probability as a function of temper- 
ature, where the MWPMA is used for decoding. The simula- 
tions are performed on the RHG lattices of the size N x N x N 
for N = 6,8, 10, 12, 14. 



The probability that the error chains belong to the 
same homology class as that of C + Vi is given by 



-"C'+V 



08) 



= exp[-P[F c . +£i (p) - Ft* 08)]], 

where f3F tot ((3) = — In -Zc'+v, (/?)] i s independent on 
the homology class i. This indicates that the homology 
class i m ; n that has the minimum free energy is the most 
likely to occur. Thus if we perform the recovery operation 
according to the estimated error chain C'+Vi min , the error 
correction succeeds with a high probability. The logical 
error probability pl is the probability that the actual 
and estimated error chains belong to different homology 
classes, namely, = 1 — Pi min - In the Higgs (ordered) 
phase, the logical error probability pl decreases exponen- 
tially with the system size. 

The free energy can be calculated approximately by 
numerical simulations, e.g., thermodynamic integration 
by use of data obtained in the Monte-Carlo simulations, 
and Population annealing (3l| with the aid of the Jarzyn- 
ski equality [32| . For the decoding method using the free 
energy to be efficient, the simulation time has to be poly- 
nomial in the system size. Although the relaxation time 
to equilibrium state takes exponentially long time in spin- 
glass phases due to a large number of metastable configu- 
rations, such spin-glass behaviors are not expected in the 
Higgs (ordered) phases. Since we perform topologically 
protected MBQC in the Higgs phase appearing below the 
critical temperature, it is expected that the error correc- 
tion using the free energy would be done efficiently [l4| . 
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